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ABSTRACT 

The Edgeworth-Kuiper belt (EKB) and its presumed dusty debris is a natural reference for extrsolar debris disks. We re-analyze 
the current database of known transneptunian objects (TNOs) and employ a new algorithm to eliminate the inclination and the 
distance selection effects in the known TNO populations to derive expected parameters of the "true" EKB. Its estimated mass is 
^EKB = 0. 12M5, which is by a factor of ~ 15 larger than the mass of the EKB objects detected so far. About a half of the total EKB 
mass is in classical and resonant objects and another half is in scattered ones. Treating the debiased populations of EKB objects as dust 
parent bodies, we then "generate" their dust disk with our collisional code. Apart from accurate handling of destructive and cratering 
collisions and direct radiation pressure, we include the Poynting-Robertson (P-R) drag. The latter is known to be unimportant for 
debris disks around other stars detected so far, but cannot be ignored for the EKB dust disk because of its much lower optical depth. 
We find the radial profile of the normal optical depth to peak at the inner edge of the classical belt, a 40 AU. Outside the classical EKB, 
it approximately follows r oc r^^ which is roughly intermediate between the slope predicted analytically for collision-dominated (r"'-^) 
and transport-dominated (r disks. The size distribution of dust is less affected by the P-R effect. The cross section-dominating 
grain size still lies just above the blowout size (~ 1 ... 2 [xm), as it would if the P-R effect was ignored. However, if the EKB were by 
one order of magnitude less massive, its dust disk would have distinctly different properties. The optical depth profile would fall off 
as T oc and the cross section-dominating grain size would shift from ~ 1 ... 2 \.im to ~ 100 \im. These properties are seen if dust is 
assumed to be generated only by known TNOs without applying the debiasing algorithm. An upper limit of the in-plane optical depth 
of the EKB dust set by our model is r = 2 x lO""" outside 30 AU. If the solar system were observed from outside, the thermal emission 
flux from the EKB dust would be about two orders of magnitude lower than for solar-type stars with the brightest known infrared 
excesses observed from the same distance. Herschel and other new-generation facilities should reveal extrasolar debris disks nearly 
as tenuous as the EKB disk. We estimate that the Herschel/PACS instrument should be able to detect disks at a ~ 1 . . . 2Mekb level. 

Key words. Kuiper belt: general - methods: statistical - methods: numerical planetary systems - circumstellar matter - infrared: 
planetary systems. 



1. Introduction 

Debris disks, now known to be ubiquitous around main- 
sequence stars, are the natural aftermath of the evolution of 
dense protoplanetary disks that may or may not result in for- 
mation of planets (see, e.g., Wvatt ,2008t lKrivovl l2010. for re- 

[ cent reviews). They are composed of left-over planetesimals and 

. smaller debris produced in mutual collisions, and it is the tiniest, 
dust-sized collisional fragments that are evident in observations 

[ through thermal radiation and scattered stellar light. 

Like planetary systems of other stars, our solar system con- 

' tains planetesimals that have survived planetary formation. The 
main asteroid belt between two groups of planets, terrestrial and 
giant ones, comprises planetesimals that failed to grow to plan- 

' ets b ecause of the strong p erturbations by the nearby Jupiter 
(e.g.. ISafronovl [19691 lwSherill .1980 ). The Edgeworth-Kuiper 
Belt (EKB) exterior to the Neptune orbit is built up by plan- 
etesimals that did not form planets because the density of the 
outer solar nebula wa s too low (e.g., Safronov 1969; Lissauer 
ll987LlKenvon & Broml ev 2008). Both the asteroid belt and the 
Kuiper belt are heavily structured dynamically, predominantly 
by Jupiter and Neptune respectively. They include non-resonant 
and resonant families, as well as various objects in transient 
orbits ranging from detached and scattered Kuiper-belt objects 
through Centaurs to Sun-grazers. Short-period comets, another 
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tangible population of small bodies in the inner solar system, 
must be gene tically related to the Kuiper belt that acts as their 
reservoir (Ou inn et al.l ll990). Asteroids and short-period comets 
together are sources of interplanetary dust, observed in the plan- 
etary region, although thei r relative contribu tion to the dust pro- 
duction remain uncertain dGriin et al.ll200lh . And this complex 
system structure was likely quite different in the past. It is ar- 
gued that the giant planets and the Kuiper belt have originally 
formed in a more compact configuration (the "Nice model", 
iGomes et al.l l2005). and that it went through a short-lasting pe- 
riod of dynamical instability, likely explaining the geologically 
recorded event of the Late Heavy Bombardment (LHB). 

As the amount of material and spatial dimensions of the EKB 
surpass by far those of the asteroid belt and the population of 
short-period comets, it is the EKB and its presumed collisional 
debris that should be referred to as the debris disk of the solar 
system. Ironically, the observational status of the solar system's 
debris disk is the opposite of that of the debris disks around other 
stars. In the latter case, as mentioned above, it is dust that can be 
observed. In the former case, we can observe the p lanetesimals, 
but there is no certain detec tion of their dust so far (iGurnett et al.l 
ll997HLandgraf et al.ll2002h . 

An obvious difference between the debris disks detected so 
far around other stars and our solar system's debris disk is the to- 
tal mass (and thus, also the amount of dust). Miiller et al. (2010l) 
for example infer several Earth masses as the total mass of the 
Vega debris disk, whereas the Kuiper belt mass is reported to 
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be below one-tenth of t he Earth mass (iBernstein et alJ l2004t 
iFuentes & Hobiianll2008h . As a resuh, were the solar system ob- 
served from afar, its debris dust would be far below the detec- 
tion limits. However, a number of debris disks around Sun-like 
stars that are coeval w ith or even older than the Sun have been 
detected. iBooth etal] (^009) analyze "dusty consequences" of a 
major depletion of the planetesimal populations in the solar sys- 
tem at the LHB phase. They point out that the pre-LHB debris 
disk of the Sun would be among the brightest debris disks around 
solar-type stars currently observed. Future, more sensitive obser- 
vations (for instance, with Herschel Space Observatory) should 
detect lower-mass disks, bridging the gap between dusty debris 
disks around other stars and tenous debris disks in the present- 
day solar system. 

Given the low mass of the dust parent bodies in the 
EKB, our debris disk is thought to fall into the category o f 
the so-called transport-dominated disks jKrivov et al.l 2000l) . 
where radial transport timescales for du st (due to the Poyn ting- 
Robertson effect, henceforth P-R, e.g. iBurns et aP Il979h are 
shorter than their coUisional lifetime. This is opposite to ex- 
trasolar debris disks that are collision-dominated (Wvatt 2005). 
Whereas the latter have b een extensively modeled b oth an- 
alytic ally and numerically ('The bault et a l.l 12003; Kriv ov et al.l 
2006t [Str ubbe & Chiang 2006; Thebauh & Augereau HooS 
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effort has to be invested into modeling of transport-dominated 
disks. The more so, as new-generation facilities like Herschel 
Space Observatory will soon open a phase when low-density, 
transport-dominated disks can be studied observationally. 

The main goal of this paper is to develop a more realistic 
model of such a tenuous dis k, exemplifie d by the EKB dust disk, 
than was done before (e.g. ISternlll995l) . We take an advantage 
that — unlike with other debris disks and unlike at the time when 
the first collisional models of the EKB dust were devised — 
more than a thousand EKB objects, acting as dust parent bod- 
ies of the solar system's debris disk, are now known. This task is 
accomplished in two steps: 

I. First, the currently known objects in the EKB are analyzed. 
In Sect. 2, we work out an algorithm to correct their distribu- 
tions for observational selection effects and try to reconstruct 
the properties of the expected "true" EKB. 
II. Second, we treat the objects of the "true" EKB as dust parent 
bodies. In Sect. 3, we make simulations of dust production 
and evolution with a statistical code, fully taking into ac- 
count collisions and the P-R effect, and present the expected 
radial and size distribution of the presumed EKB dust. In 
Sect. 4, we model the spectral energy distribution (SED) of 
the simulated EKB dust disk and compare it to the SEDs of 
other debris disks. We also compare these results with the 
detection limits of the Herschel/PACS instrument. 

Our results are summarized in Sect. 5 and discussed in Sect. 6. 



2. Planetesimals in the Kuiper belt 

2.1. Observations and their biases 

The EKB was predicted more than fifty years ago by Edgeworth 
and Kuiper and i t took fourty year s until the first object, QB 1, 
was discovered jJewitt et al.lll992l) . More than 1260 transnep- 
tunian objects (TNOs) orbiting the Sun beyond the orbit of 
Neptune have been discovered. Table [T]lists most of the surveys 
published so far, in which new TNOs have been discovered, and 
key parameters of these surveys. One parameter is the area Q. on 
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Jewitt&Luu(1995) 
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Gladmanet al. (1998) 
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Gladmanet al. (1998) 
Jewitt etal. (1998) 
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Truiillo et al. (2000) 
Al len etal. (20011 * 
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Truiillo et al. (2001b)* 
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Elliot et al. (2005)* 
Larsen et al. (2007)* 
Eraser et al. (2008) 
Euentes & Holman (2008) 

Euentes et al. (2009) 
Eraser & Kavelaars (2009) 



Table 1. A list of campaigns where TNOs were found. The sky 
area covered (O), the number of the objects discovered (A^), the 
limiting magnitude (niso), an estimated half-opening angle e, 
and the ecliptic longitude coverage a are given. Papers that pro- 
vide enough data for objects discovered in that survey to identify 
them in the MPC database are marked with an asterisk. 



the sky searched for TNOs. Another one is the limiting magni- 
tude ;«5() that corresponds to the detection probability of 50%. As 
the detection probability drops rapidly from 100% to zero when 
the apparent magnitude m "crosses" m^o, we simply assume that 
an object will be detected with certainty if m < m^Q and missed 
otherwise. Finally, the maximum ecliptic latitude s covered by 
each survey is listed. Where it was not given explicitly in the 
original papers, we estimated it to be Vn/2, assuming that the 
surveyed area was centered on the ecliptic. Table[T]shows that all 
campaigns can be roughly divided into two groups: deeper ones 
with a small sky area covered ("pencil-beam" surveys) and shal- 
lower ones with a larger area, but a smaller limiting magnitude. 

The orbits of TNOs are commonly characterized by six or- 
bital elements: semi-major axis a (or perihelion distance q), ec- 
centricity e, inclination /, argument of pericenter co, longitude of 
the ascending node Q, and mean anomaly M. In addition, each 
object itself is characterized by the absolute magnitude H, which 
is defined as the apparent magnitude the object would have if it 
was 1 AU away from the Sun and the Earth, and depends on the 
object radius and albedo. We take the orbital elements and the 
absolute magnitudes of all known objects from the Minor Planet 
Center (MPC) databas43 rather than from discovery papers listed 
in Tab.[Tl because the MPC data include follow-up observations 
and thus provide a better precision. 

The planet formation theory implies that the TNO orbits 
strongly concentrate towards the ecliptic plane. Accordingly, 
in order to increase the detection probability, the majority of 
the observations were made near the ecliptic, and only a few 



' http://www.cfa.harvard.edu/iauJlists/TNOs.html [Last accessed 12 
October 2009] 
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Fig. 1. Known TNOs in the a-e plane (top) and a-i plane (bottom). Different groups are shown with different symbols: 865 classical 
objects with dots, 235 resonant TNOs with pluses, and the remaining 160 scattered objects with crosses. Solid lines separate classical 
and scattered objects in our classification. One object with / = 103.4° is outside the lowest panel, but is marked with arrows in both 
panels. Notice that the linear scale turns into a logarithmic scale at 100 AU. 



surveys covered high ecliptic latitudes. Our sample, given in 
Tab. [1] contains surveys with e up to 10° ( e.g. T r ujillo et al, 
200 laF iT ruiillo & Brown 2003; lEUiot et all lIooBi: iPetit et alj 



2006t iLarsen et aL.2007.) . How ever, TNOs w ith higher orbital 



inclinations exist as well. Since iBrownl (1200 lb it is known that 
the inclination distribution of TNOs has a second maximum at 
higher inclinations. Several objects with very high inclinations, 
including one with / - 103.4°, were detected. Clearly, the fact 
that observations are done near the ecliptic plane decreases the 
probability to detect such objects, because it is only possible 
twice per orbital period, close to the nodes. Thus there is an ob- 
vious selection effect in favor of TNOs in low-inclination orbits 
that needs to be taken into account. Equally obvious is another 
selection effect, which is that objects are predominantly discov- 
ered at smaller heliocentric distances. This reduces the proba- 
bility to discover TNOs with large semi-major axes and high 
eccentricities, because these are too faint all the time except for 
the short period of time when they are near perihelion. 



2.2. Classification of TNOs 

Many classifications of TNOs into "classical", "resonant", "ex- 
cited", "scattered", "detached" etc. groups have been proposed, 
based on the orbital elements and taking into account dynam- 
ical arguments (e. g., Jewitt et al. 1998; Chiang & Brown 1999; 
lLuu& Jewit3l2002i: .Delsanti & Jewitt. ,2006i: Jewitt et al.li2009L 



among others). Classifications by different authors are similar, 
but not identical. In this paper we use the following working 
classification: 

1. Resonant objects (RES): objects in a mean-motion com- 
mensurability with Neptune, where we only consider the 
three most prominent first-order resonances 4:3, 3:2 and 2:1. 
To identify th e objects as resonant, we use the resonance 
"widths" from iMurrav & DermottI (l2000h . For example, the 
width of the 3 : 2-resonance at e = 0.1 is Aa = 0.012 AU. 
The width increases with increasing eccentricity and with de- 
creasing distance to Neptune. 

2. Classical Kuiper Belt (CKB) objects: objects with a < 
50 AU, which are neither resonant nor Neptune-crossers (q > 

'JNeptune)- 

3. Scattered disk objects (SDO): objects with a > 50 AU, as 
well as Neptune-crossers (a > flNeptune and q < flNeptune)- 

Figure [T] depicts all known TNOs, using different symbols for 
each of the three groups. This classification is intentionally made 
simpler than many others in common use, in order to facilitate 
the analysis below. For instance, in our classification, the known 
detached objects (those with perihelia well outside Neptune's or- 
bit, q > 40 AU) fall into the "scattered" category. On any ac- 
count, the parameters of the entire EKB and its dust that we will 
derive will not depend on the way in which the TNOs are clas- 
sified into various groups. On the other hand, this classification 
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roughly reflects difi'erent formation history of different popula- 
tions in the EKB, as well as diff'erent modes of their gravitational 
interaction with Neptune at present. 

2.3. Debiasing procedure 

Because of the obvious selection efifects of inclined and faint 
objects, statistical models were developed to estimate a t rue dis- 
tributi on of orbital elements and numbers of the TNOs. iBrownl 
1200 11) calculated an inclination distribution. He assumed circu- 
lar orbits and derived a relation between the inclination and the 
fraction of an object's orbit that it spends at low ecliptic latitudes. 
iDonnisoiil (|2006) calculated the magnitude distribution for the 
classical, resonant, and scattered objects for absolute magnitudes 
H <1, using maximum likelihood estimations, and showed that 
the samples were statistically different. 

Here we propose another debiasing method to estimate the 
"true" distribution of the TNOs, based on the obervational sur- 
veys listed in Table [1] We start with calculating the probability 
to find an object with the given parameters {a, e, i, co, Q, H] for 
each given survey. To this end, we estimate the time fraction of 
the object's orbit that lies within the maximum ecliptic latitude s 
covered by the survey, as well as the fraction of the orbit which 
is observable at the given limiting magnitude mso, and find the 
intersection of the two orbital arcs. Once the probability to de- 
tect the object in each of the surveys has been calculated, we 
compute the probability that it would be detected at least in one 
of the surveys made so far Finally, we augment the number of 
objects with that same orbital elements as the object considered 
to a 100% probability, e.g., an object with 20% probability is 
counted five times. 

We now explain this procedure in more detail. The first effect 
is the "inclination bias". In calculating the orbital arc that lies 
in the observable latitudinal zone, we make the assumption that 
we observe from the sun. The observable area on the sky is thus 
a belt \b\ < e, where b is the heliocentric ecliptic latitude. The 
orbit crosses the boundary of the observed belt, \b\ - s, at four 
points. At these intersection points, the true anomaly (p takes the 
values 
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sin^ / 



- O). 



(1) 



Due to our approximation that we observe from the Sun, the lon- 
gitude of the ascending node does not appear in this formula. 

The second effect is the "distance bias" (or "eccentricity 
bias"). The maximum distance at which an object is detectable 
is given by Clrwin et al..,1995 ) 
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Then we combine both observability constraints, from incli- 
nation (Eq. [1]) and eccentricity (Eq. |2|, into one, to find the or- 
bital arc (or arcs) that lie both in the observable latitudinal belt 
and the observable sphere. Typical geometries are sketched in 
Fig. 121 assuming that the pericenter is inside the observing latitu- 
dinal belt. The intersection points of the orbit with the visibility 
sphere r = r^ax are denoted by E^, those with the visibility circle 
\b\ = ±s by Ik (indices k increase with increasing true anomaly). 
The point Ei can lie before /i, between Ii and I2, or between 
I2 and I2, giving three possibilities. On the other hand, the point 
E2 can reside between I2 and 73, between I3 and I4, or after I4. 
This yields 3x3 = 9 possibilities in total, denoted by I... IX. 
Remember that the position of the observable arcs in Fig.|2]cor- 
responds to the case where the perihelion is inside the observable 
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Fig. 2. Observable arc(s) of a TNO orbit that satisfy the dis- 
tance and the inclination restrictions. Smaller points denoted 
with £1 and E2 are intersection points of the orbit with the sphere 
r - /"max- Bigger points I\ . . .I4 are intersection points of the or- 
bit with the circles b - +e on the sky. Solid and dashed arcs 
represent observable and non-observable parts of the orbit, re- 
spectively. 



latitudinal belt. If it is outside, the observable arcs change, giv- 
ing rise to another set of nine possibilities. Thus there are 18 
possibilities in total. Additionally, there are special cases. One is 
/ < e, where the entire orbit is inside the observable belt, so that 
the points /i ... 74 do not exist. Others are where the entire orbit 
is inside or outside the observable sphere, so that the points Ei 
and £2 do not exist. 

As an example, we take ellipse number III. The object starts 
at the pericenter (where it is visible) and moves toward I\. 
Between I\ and I2, it is outside the observed latitudinal belt and 
is invisible. Although it has a sufficiently low ecliptic latitude 
up to 73, it is only detectable up to Ex, because it gets too faint 
beyond that point. Between 73 and 74 the object is too far from 
from the ecliptic, and it stays outside the limiting sphere until it 
reaches E2. Starting from E2, the object is visible again. 

Having found the observable arcs, we compute the fraction 
/ of the object's orbital period it spends in these arcs for a given 
survey. If the survey had a full (360°) coverage of the ecliptic 
longitudes, that fraction would directly give us the probability to 
detect the object. However, the survey is confined to the longi- 
tude range with a certain width a. This width is often given in 
the papers, and where it is not, we simply take a - yfO. - 2e, 
where Q. is the observed sky area. Thus the detection probabil- 
ity is / multiplied by q'/360°. This estimation assumes that all 
the observations are done within a period of time that is much 
shorter than the orbital period of the TNO, so that its proper 
motion can be neglected. This is the case for all the surveys 
we consider However, the same estimation assumes a uniform 
azimuthal distribution of the TNOs. For plutinos, for instance, 
this is no longer valid, as they concentrate preferably in two az- 
imuthal zones ahead and behind Neptune's location. Thus our 
algorithm may underestimate the detection probability of reso- 
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nant objects, at least in surveys targeted at parts of the ecliptic 
where such objects are more numerous. 

In this way, for each of the known TNOs, we can calculate 
the detection probability in any survey. We then calculate the 
probability P,- that an object / would be detected in any of the 
A^surveys survcys: 

"surveys 

Pi^i- n n-Pik), (3) 

k=l 

where Pik is the probability to detect an object / in a survey k. 
The advantage of equation (3) is that 1 - f ^ gives the probability 
not to detect an object, "shallower" surveys make little contribu- 
tion to the product and thus to the total detection probability of 
very faint objects. Therefore, it is deep surveys that dominate the 
result for faint objects. 

Given the discovery probability Pi of a given object, we can 
augment the observed Kuiper belt to the "true" one by counting 
that object Pr' times. In other words, we debias the observed 
Kuiper belt by setting the number of TNOs with the same orbital 
elements as the known object to Pr'. 

The number of surveys in Tab. [T] is A^surveys - 23. However, 
only nearly half of the 1260 TNOs contained in the MPC 
database were found in these campaigns. Another » 600 objects 
were discovered in other observations, some serendipitously in 
surveys that did not aim to search for TNOs. The circumstances 
of those observations have not been been published in all cases. 
What is more, even for the campaigns listed in Tab. [1] it is prob- 
lematic to identify which particular set of a; 600 objects out 
of 1260 in total was found in those surveys. Indeed, the papers 
that give a specific, identifiable list of newly discovered objects 
(marked with an asterisk in Tab. [U only cover » 400 TNOs. 
We do not know under which circumstances the remaining two- 
thirds of the TNOs were detected. In other words, there is no 
guarantee that the parameters of those unknown surveys (wgo, e 
etc.) are similar to those listed in Tab. [1] Furthermore, some of 
the surveys in our list may not have reported their discoveries to 
the MPC. As a result, it is diflicult to judge how complete the 
MPC database is. We can even suspect that there have been sur- 
veys not listed in Tab. [T] that have discovered TNOs not listed 
in the MPC. Therefore, it does not appear possible to compile 
a complete version of Tab. [T] that would cover all known TNOs 
and all discovery observations (together with their Q, mso, and 
e). Nor is it possible to get a complete list of all known TNOs 
together with their orbital elements, along with information in 
which particular survey each of the known TNOs was discov- 
ered. To cope with these difficulties, we make two assumptions. 
First, we assume that the surveys listed in Tab. [1] are represen- 
tative of all surveys that discovered TNOs. Second, we assume 
that, conversely, the TNOs listed is the MPC are respresentative 
of all the objects discovered in surveys listed in Tab. [T] These 
two assumptions represent the main shortcoming of our debias- 
ing approach. 

To check them at least partly and proceed with the debias- 
ing, we employed two dififerent methods. In the first method, we 
have randomly chosen 600 TNOs out of the full list of known 
objects and assumed that it is these objects that were discovered 
in the campaigns listed in Tab. [T] We tried this several times for 
different sets of 600 TNOs and found that the results (e.g., the el- 
emental distributions and the total mass of the "debiased EKB") 
are in close agreement. In the second method, we have made an 
assumption that another set of 23 similar surveys with similar 
detection success rate would have likely led to a discovery of 
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Fig. 3. The detection probability of all known TNOs (as func- 
tion of their orbital inclinations) in a fiducial survey with a full 
coverage of a belt on the sky within e = 5° ecliptic latitude and 
the limiting magnitude of m^Q - 25 mag. The curve is the formal 
detection probability of objects in circular orbits in the mso - oo 
limit, but it approximates well the detection probability of many 
known TNOs in eccentric orbits in our fiducial survey. Objects 
which are below the curve are either those affected by the dis- 
tance bias or have arguments of pericenters which are outside 
of our viewing field. Objects above the curve correspond to rare 
cases where the orbital eccentricity is high, aphelion lies in the 
observable belt, and the object is not too faint even near the aphe- 
lion. These are mostly scattered objects. 



all known TNOs. So we simply counted each survey twice and 
replaced Eq. Q by 

N 

' ' surveys 

A' = 1 - n (1 - 'P''^)' ■ (4) 

k=\ 

Again, the results turned out to be very close to those found with 
the first method. 

Figure [3] illustrates the probabilities Pjk to observe known 
TNOs in a fiducial survey with mgo = 25 mag, a latitudinal cov- 
erage of e = 5°, and a longitudinal coverage of 360°. Let us start 
with an artificial case where all objects are in circular orbits. If 
they were bright enough to be observed (or equivalently, in the 
limiting case niso - oo), they would all lie on the curve overplot- 
ted in Fig. [3] In particular, their detection probability would be 
100% for / < e, and it would be e/90° = 5.6% for / = 90°. If 
they are too faint for detection, their detection probability will 
be zero regardless of the inclination. The case of eccentric orbits 
is more complicated. Then, the vast majority objects still are on 
the curve but, as seen in the figure, there are many that lie be- 
low. Either these are objects whose pericenter is outside of the 
latitudinal belt \b\ < e or these are objects that cannot be ob- 
served over the entire orbits, even when they have sufficiently 
low ecliptic latitude, because in some low-latitude parts of their 
orbits they are too faint to be visible. In fact, a mixture of both 
cases is typical. Finally, a few objects lie above the curve. These 
are rare cases of objects in highly-eccentric orbits, whose aphe- 
lia fall into the observable latitudinal belt, and whose apocentric 
distances are not too large. Such objects spend much of their or- 
bital period near aphelia and are detectable there, which raises 
their detection probability. 
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Fig. 4. The final detection probabihty of the known TNOs, cal- 
culated with Eq. Included are all surveys from Table[T] 



Although the average detection probability in Fig. [3] is quite 
high, this only holds for a complete coverage of the \b\ < s 
band on the sky. In reality, only a limited range of the eclip- 
tic longitude is covered. The resulting detection probability P, 
of all known TNOs in all 23 surveys, calculated with Eq. ^ 
that takes into account actual latitudinal coverage of the obser- 
vational campaigns, is plotted in Fig.|4] Typical values are within 
~ 20% for near-ecliptic orbits and drop to a few percent for in- 
clinations above 10°. 



TNOs are also large, see middle panels in the bottom rows of 
Figs. |5] and |6]l. These conclusions should be taken with caution, 
because the statistics of scattered objects is scarce and their de- 
biasing factors are the largest. 

The eccentricity distribution in Figs. |5] and |6] shows moder- 
ate values (e < 0.2) for the classical belt and reveals a broad 
maximum at e ^ 0.1 . . .0.3 for the resonant objects. The max- 
imum for the scattered objects appears to be located around 
e ~ 0.5... 0.6. 

As far as the inclination distribution (right panels in Figs J5] 
and|6) l is concerned, our analysis confirms the result by Brow^ 
('2001) who indentified two distinct subpopulations in the clas- 
sical Kuiper belt, a cold one with low inclinations and a hot 
one with more inclined orbits. The maxima of 0° ... 5° and 
20°... 25° that we found are consistent with his results of 
2°6!5;:^I and 17° ± 3°. 

The inclination distribution of the resonant o bjects reveals a 
broad maximum around a: 15°. For comparison, iBrownl (1200 Ih 
found a maximum at 10°2^[ g„. A second maximum visible at 
/ - 30° . . . 35° in the number distribution (Fig.|5]l is due to small 
objects with a large debiasing factor, which are still big enough 
not to fall under the i < 75 km criterion. That is why in Fig. |6] 
the same peak is barely seen. 

A clear difference between the number and mass distribu- 
tions can be seen in the bottom right panels of the two fig- 
ures, too, which show the inclination distribution of the scat- 
tered objects. A large number of scattered TNOs can be found 
at 25° . . . 30° (Fig. 0, whereas their mass peaks at 15° 20° 
(Fig. |6]l. Interestingly, a recent paper by iGulbis et al.l (1201 Oh 
yielded 19°1^3^„, which is close to the maxima we find here. 



2.4. Orbital element distributions of the Kuiper belt objects 

Having applied the debiasing procedure, we compared and ana- 
lyzed the distributions of orbital elements of the known and the 
"true" EKB — separately for each class. 

Figures |5] and |6] show the distributions in terms of num- 
bers (for s > 75 km) and masses (for s < 400 km) of objects 
per element's bin before and after debiasing. The distribution in 
terms of numbers depicted in Fig. |5] emphasizes smaller, more 
numerous, TNOs. It is directly related to observational counts 
of TNOs and is also useful to alleviate comparison with sim- 
ilar work by the others. In contrast, the distribution of TNO's 
mass in Fig. |6]is dominated by larger objects. It demonstrates 
more clearly where the wealth of the EKB material is located, 
which aids placing the EKB in context of extrasolar debris disks. 
Objects with s < 75 km were excluded from Fig. |5] because 
detections of the smallest objects are the least complete, which 
would lead to a highly uncertain, distorted distribution. And con- 
versely, we excluded the biggest objects with s > 400 km from 
Fig|6]to avoid large bin-to-bin variations stemming from a few 
individual rogues. Failure to do this would lead, for example, 
to a pronounced peak in the eccentricity distribution of resonant 
objects at e = 0.15 . . . 0.20 produced by a single object, Pluto. 

As seen in Figs. |5] and |6] for the classical Kuiper belt, de- 
biasing increases the total number and mass of objects, but the 
position of the maximum remains at a » 44 AU. The same holds 
for the semimajor axis distribution of resonant objects, whose 
peaks are preserved at known resonant locations. In contrast, for 
the scattered objects, here are indications that a substantial unbi- 
ased population with larger semimajor axes of 80 AU ... 120 AU 
might exist. Some of them may be "detached" (q > flNeptune), 
while some others may not (since the eccentricities of these 



2.5. Albedos and sizes of the Kuiper belt objects 

To e stimate the TNO sizes , we employed the V-band formula 
from lKavelaars et al.1 (120091) : 

H = mo + 42.38 - 2.5 lg(4ps^), (5) 

where = -26.74 mag is the apparent V-magnitude of the Sun, 
p the albedo and s the radius of an object in kilometers. Solving 
for radius, we find 

i = 671.5 X km. (6) 

With th is equation and albedo measurements from [ Noll et al.l 
(12004 : IStansberrv et"an (l2008l) : iBrucker et al.l ilOO^L we cal- 
culated the radius of objects with known albedo (Fig.|7]i. 

Albedos inferred for a handful of big objects with H ^ 3 
turned out to be high, which is indicative of a strongly reflecting 
surface material. For instance, the surface of Haumea wa s found 
to be c overed with > 92% pure water ice (Pinilla-Alons o et aH 
l2009h . iDumas et all (l2007l) reported for Eris 50% methane ice 
on its surface along with nitrogen and water ices, and ice tholin. 
Smaller objects are coated with darker carbonaceous layers, so 
their albedo is lower. Note that objects between 6 < H <1 have 
a very strong scatter, the reason for that being unknown. Albedos 
of the smallest TNOs with 1 < H < 9 are typically close to 
~ 0.05, and there have been no measurements beyond H - 9. 
However, since the EKB is known to act as a reservoir of short- 
period comets, we can use the measurements of cometary nuclei 
with sizes of ~ 1-10 km as a proxy for the reflectance properties 
of the smallest TNOs. (The obvious caveat is that comets may 
have altered their original surface properties as a result of their 



Vitense et al.: The Edgeworth-Kuiper debris disk 



7 




20x observed 



35 40 45 50 0.2 0.4 0.6 0.8 1 10 20 30 40 




CO 

c 

O 
W 
CD 




M 



35 40 45 50 0.2 0.4 0.6 0.8 1 10 20 30 40 




I expected ^^^H 
50x observed 
; J 




40 60 80 100 120 140 160 180 200 0.2 0.4 0.6 0.8 1 10 20 30 40 50 

a[AU] e i(1 

Fig. 5. Distribution of classical (top row), resonant (middle row), and scattered objects (bottom row), in terms of numbers of objects. 
Left column: semimajor axes, middle: eccentricities, right: inclinations. Black and grey bars in each panel represent the expected 
(unbiased) and observed populations, respectively. The numbers of the observed TNOs are magnified by 20 (classical and resonant 
objects) and 50 (scattered objects) for better visibility. Numbers are given in 1000 for intervals with a width of Aa = 1 AU (classical 
and resonant), Aa - 10 AU (scattered objects), Ae = 0.1 and Ai - 5° for all populations. 



long residence in the inner solar system.) The typical a^lbedo val- 
ues of the nuclei range from 0.02 to 0.06 dLamv et al.ll2004 . 

On these grounds, to eliminate the dependence on albedo 
(which is not known for most of the Kuiper belt objects) from 
Eq. (|6]l, we have fitted the sizes of the TNOs with known albedo 
by an exponential function at H < 6 and assumed p = 0.05 for 
all TNOs with H > 6. This yielded a formula where s is only a 
function of H: 



s = 926 X 10 



-0.1 19H 



km {H <6) 



and 



3000 X 10"°-^^ km 



(H > 6). 



(7) 



(8) 



The smallest object found so far is a scattered object with H - 
15 mag which corresponds to a size of only ^ = 3 km. The small- 
est resonant object has a radius of s - 9.9 km (H = 12.4 mag) 
and the smallest classical one has s = 12.5 km. 



2.6. Mass of the Kuiper belt 

To translate the TNO sizes into mass requires assumptions 
for the bulk density. We took the commonly used value of 
p = 1 gcm"^. This assumption is in accord with the values of 
(()^6^. .2.3)gcm"^ found for a few individual TNOs (Lacerd^ 
l2009l) . The resulting mass and number of objects in resonant, 
classical, and scattered populations and in the entire Kuiper 
belt are listed in Tab. [2l The deduced "true" ma sses are sev- 
eral times higher than in lFuentes & Ho lman 72008^ who inferred 
McKB = (0.008 + 0.001)Me, MsDO = 0.010!°;°2iMe, with a to- 
tal of Mtot = 0.020^Q QQ3Me. However, they considered the mass 
within ±3° around the ecliptic. Since we investigated the full 
range of ecliptic latitudes, we deem the results consistent with 
each other. 

One issue about the deduced mass of the entire EKB and 
its populations is the influence of the uncertainties of the or- 
bital elements inferred from the observations. In many cases, 
the elements are known only roughly, and some of them are not 
known at all. For 134 out of 865 known classical objects, for 
instance, the observed arc was so short that a circular orbit was 



8 



Vi tense et al.: The Edgeworth-Kuiper debris disk 




Fig. 6. Same as Fig.|5] but in terms of mass contained in TNOs. 
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Fig. 7. Absolute magnitude - radius relation for objects with 
known albedo. The biggest objects are labeled with their names. 
The thick solid line is a fit to this relation, Eqs. O-®. Thin lines 
correspond to equal albedos of p = 0.05; 0.10; 0.20; 0.50; 1.00. 



assumed by the observers (these objects are clearly visible as an 





before debiasing 


after debiasing 


McKB+RES [O.OOlMe] 


4.4 


52 


MsDO [O.OOlMe] 


2.2 


65 


M,„„i [O.OOlMe] 


6.6 


117 


A^CKB+RES {s < 75 km) 


1100 


39700 


A'sDO {s < 75 km) 


160 


45200 


Motai (s < 75 km) 


1260 


85000 



Table 2. Masses and numbers of objects in the Kuiper belt. CKB 
is the classical Kuiper belt, RES are the resonant objects, and 
SDO denotes the scattered disk objects. 



e - stripe in the upper panel of Fig. [T]). How could a change 
in the orbital elements of an object affect the debiasing proce- 
dure and the final estimates of the parameters of the "true" EKB? 
Obviously, if a true value of one of the three elements of a TNO 
(a, i, or the absolute magnitude H) is larger than the one given 
in the database, the detection probability will be overestimated 
and the estimated number of similar objects in the "true" EKB 
underestimated. The eccentricity plays a special role in this case. 
Increasing it would not automatically lead to a lower detection 
probability, the pericenter distance decreases while the apocen- 
ter increases, so that the total detection probability depends also 
on a. However, a combined variation of two or more elements 
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o-[%] 


AfcKB 


Mres 


A^CKB+RES 


5 


46.1 ± 1.6 


8.7 ±0.6 


54.8 ± 1.7 


10 


48.3 + 2.2 


9.1 ±0.7 


57.3 ± 2.3 


15 


53.3 ±3.5 


9.7 ±0.8 


63.0 ± 3.6 



Table 3. Masses of objects in the Kuiper belt, as a function of 
the assumed relative accuracy cr, with which orbital elements of 
TNOs were deduced from observations. Abbreviations and units 
are as in Tab. |2] 



may alter the results in either direction. As an example, let us 
consider a scattered object with a - 1057 AU and e - 0.977, 
which has a pericenter distance of q - 24.3 AU. Decreasing, 
for instance, both a and e by 5% would lead to a pericenter at 
q -12 AU, which would result in a significantly lower detection 
probability and therefore in a higher contribution of that object to 
the estimated total mass. In contrast, we may consider a classical 
object with a - 40 AU and e - 0.2, which cannot be observed 
near the apocenter. Again, decreasing both values by 5% would 
now reduce the aphelion distance to detectable values, so that the 
detection probability would increase. 

From published observational results, we assume 5-10% 
as a typical error for the orbital elements. To quantify possi- 
ble effects, we used the following Monte-Carlo procedure. We 
assumed that the orbital elements and the absolute magnitude 
[a, e, i, H] are known with a certain relative accuracy cr (for sim- 
plicity, the same for all four elements). Then, we randomly gen- 
erated [a, e, i, H] -sets for each of the known TNOs assuming 
that each element of each object is uniformly distributed be- 
tween X- cr and x + cr, where x is the cataloged value. For this 
hypothetical EKB, the debiasing procedure was applied and the 
expected masses and number of objects in the "true" Kuiper belt 
were evaluated. This procedure was repeated 10,000 times (for 
10, 000 realizations of the observed Kuiper belt, that is to say). 
From these calculations we excluded the scattered objects, be- 
cause varying their orbital elements would lead to extremely low 
detection probabilities. The results for several cr values between 
5% and 15% are listed in Table[3] It is seen that the effect of the 
uncertainties of the orbital elements on the global parameters 
of the "true" EKB is moderate except for the SDOs. Since the 
SDOs have "extreme" orbital elements, compared to the classi- 
cal belt, an uncertainty of, e.g., 25% would alter the total mass 
by a factor of 2. 

Interestingly, the net effect of the increasing cr is that the 
mean TNO detection probabilities slightly decrease, which leads 
to somewhat higher estimates for the mass of the EKB popula- 
tions and the whole Kuiper belt. 
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Fig. 8. Cumulative numbers of the expected and known Kuiper 
belt objects with a Ai = 5 km resolution. For an easier orienta- 
tion, straight lines show the slopes that would coiTespond to the 
differential size distribution indexes of ^ = 3, 4, and 5. 

Our results are largely consistent with previous determina- 
tions (Table |4]l. For the CKB, for instance, the range between 
3.6 ± 0.1 ( Chiang & Brown 1999) and 4.8+[J^ dGIad man et akl 
was reported. In this comparison, one has to take into 
account that different authors dealt with somewhat different 
size intervals. Chiang & Brown (1999) considered objects be- 
tween ( 50 . . . 500) km, Gladman et al. (2001) and TruUll o et al.l 
(12001 ah between (50 . . . 1000)km, and Donnison (2006h be- 
tween (120. . .540) km {1 > H > 2). For the SDOs, our result s 
are also consistent within the error bars with iDonniso iil (I2006h . 
However, for the resonant objects our result departs from his ap- 
preciably. 



CKB 


RES 


SDO 


reference 


4.3 ± 0.2 


5.1 ±0.1 


2.8 ±0.1 


this paper 








Gladman et al. (199^ 


4.0 ± 0.5 






Jewitt etal. (1998) 


3.7 ± 0.2 






Luucfe Jewitt (1998) 


3.6 ±0.1 






Chiang & Brown (199J) 


4.4 ± 0.3 






Gladman et al. (200r) 








Truiillo etal. (2001a) 


4.05 ± 0.2 






Bernstein et al. (2004) 


3.97 ±0.15 


3.30 ± 0.37 


3.02 ± 0.32 


Donnison (2006) 



Table 4. Size distribution index of the Kuiper belt populations. 



2. 7. Size distribution of tlie Kuiper belt objects 

We now come to the size distribution of KBOs. The exponents 
q of the differential size distribution N(s)ds oc ds after debi- 
asing were derived with the size-magnitude relation (|7]i-([8]l. In 
doing so, we have chosen the size range 50 km < s < 170 km 
(8.9 > // > 6), and we determined the size distribution index 
separately for different populations of TNOs and their combina- 
tions. For the CICB, the result is ^ = 4.3 + 0.2. The resonant 
objects reveal a steeper slope of 5.1 + 0.1, with plutinos (in 3:2 
resonance with Neptune) having 5.3 ± 0. 1 and twotinos (2: 1 res- 
onance) having 4.0 ±0.1. This results in 4.4 + 0.2 for classical 
and all resonant objects together. In contrast, the scattered ob- 
jects have 2.8 + 0.1. Altogether, we find 3.6 + 0.1 for the entire 
EKB (classical, resonant, and scattered TNOs). 



Figure [8] shows cumulative numbers of the expected Kuiper 
belt ob jects larger than a given size. In agreement with Donnison| 
(l2006h . the profile flattens for objects s < 60 km (H < 7). The 
break in the size distribution at radii of several te ns of kilometers 
reported by so me authors , e.g., at s k 30 km bv lBernstein et al.l 
(2004) and Fra'ser ('2009') can neither be clearly identified nor 
ruled out with our debiasing algorithm. 

3. Dust in the Kuiper belt 

3.1. Simuiations 

We now move on from the observable "macroscopic" ob- 
jects in the EKB to the expected debris dust in the transnep- 
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tunian region. We employ the technique to follow the size 
and radial distribution of solids (from planetesimals down to 
dust) in rotationally-symmetric debris disks, developed in previ- 
ous papers (iK rivov et al. 2000, 2005; Kri vov et al.li2006i l2008t 
iLohne et afluOOS; .MuUer et al.. ,2010.) . Our numerical code, 

solves the Boltzmann- 



ACE (Analysis of Collisional Evolution), 
Smoluchowski kinetic equation over a grid of masses m, perias- 
tron distances q, and orbital eccentricities e of solids. It includes 
the effects of stellar gravity, direct radiation pressure, as well as 
disruptive and erosive collisions. Gravitational effects of plan- 
ets in the system are not simulated with ACE directly. Since we 
do expect Neptune to affect the dust disk in the EKB region in 
several ways, we will discuss this later in Sect. 5. The code out- 
puts, among other quantities, the size and radial distribution of 
disk solids over a broad size range from sub-micrometers to hun- 
dreds of kilometers at different time steps, and the code is fast 
enough to evolve the distribution over gigayears. 

As explained in Sect. 1, an effect of particular importance in 
transport-dominated disks is the Poynting-Robertson drag. The 
latter is now implemented in ACE through an appropriate dif- 
fusion term in the space of orbital elements, comi ng from the 
clas sical orbit-averaged equations for q and e (e.g. iBurns et al.l 
Il979i) . To clearly see the role played by the P-R drag, we ran 
ACE twice, with and without P-R. Figure |9] illustrates the dis- 
tribution of the perihelion distance q versus orbital eccentricity 
e of dust grains. This distribution was calculated both without 
(top panel) and with P-R effect (bottom panel). The location of 
the main belt is recognizable as a dark grey region in each of the 
plots. Clearly visible is the dual role played by the P-R drag. On 
the one hand, it lowers the pericentic distances of dust particles, 
filling the lower parts of the bottom panel. On the other hand, it 
circularizes the orbits of dust particles. This can be seen, for ex- 
ample, as a slight concentration of particles, whose pericenters 
are located inside the main belt, towards the e = line. 

To model the collisional evolution of the EKB, as well as the 
distribution and thermal emission of the EKB dust, we have to 
assume certain optical and mechanical properties of dust. This, 
in turn, necessitates assumptions about its chemical composi- 
tion. The surface composition of a number of b right EKB ob- 
jects has been measured (see, e.g., iBarucci eTaLI i2008. for a 
recent review); see also discussion in Sect. 2.5. These objects 
turned out to have surfaces with very different colors and spectral 
reflectances. Some objects show no diagnostic spectral bands, 
while others have spectra showing signatures of various ices 
(such as water, methane, methanol, and nitrogen). The diversity 
in the spectra suggests that these objects represent a substan- 
tial range of original bulk compositions, including ices, silicates, 
and organic solids. A single standard composition that could be 
adopted to represent "typical" EKB dust grains is therefore dif- 
ficult to find. For the sake of simplicity, for the collisional sim- 
ulation with ACZi we choose an ideal material with p - \ gcm"^ 
and geometric optics, leading to the radiation pressure efficiency 
Qpr - 1 . Next, an important property that we need for the col- 
lisional simulations is the critical specific energy 2^, which is 
the ratio of impact energy and mass of the target. It is calculated 
by th e sum of two power laws (IKrivov et al.ll2005t iLohne et"al] 
l2008l and references therein) 



(9) 



where the first and the second term represent the strength and the 
gravity regime, respectively. We took values thought to be typi- 
cal of low-temperature ice: A , - 10^ ergg" 
3b, = -0.37 and 3b, = 1.38. 




-2 



eccentricity 




',Aj = 2x10^ ergg 



eccentricity 



Fig. 9. Phase-space distribution (in e,q'-plane) of dust main- 
tained by the expected EKB. Top: without P-R effect, bottom: 
with P-R effect. Pericenter distances are in AU. Eccentricities 
e < -1 correspond to anomalous hyperbolas open outward from 
the Sun. Grayscale gives the total cross section of dust, inte- 
grated over all sizes (in arbitrary units). 



In the collisional simulation with ACE, we refrain from test- 
ing various possible material compositions of dust, for the fol- 
lowing reasons. First, each ACE run requires up to a few weeks 
of computing time in parallel mode on 8-16 kernels. Second, we 
would need to consistently modify both the optical constants of 
the assumed material (that determine the strength of the radi- 
ation pressure force through Qp, ) and its mechanical properties 
that control the collisional cascade. The latter would be problem- 
atic because of the lack of experimental data, for instance on 2^, 
for non-icy materials. Nevertheless, in Sect. 4, we will test the 
influence of various materials on the resulting termal emission 
of the EKB dust to get a rough idea of the material dependence 
of the simulation results. 

In the ACE simulations, we used the following size- 
pericentric distance-eccentricity mesh. The minimum grain ra- 
dius was set to 0.1 (im and the variable mass ratio in the adjacent 
bins between 4 (for largest TNOs) and 2.1 (for dust sizes). The 
pericenter distance grid covered 41 logarithmically-spaced val- 
ues from 4 AU to 200 AU. The eccentricity grid contained 50 
linearly-spaced values between -5.0 and 5.0 (eccentricities are 
negative in the case of smallest grains with /3 > I, whose orbits 
are anomalous hyperbolas, open outward from the star). The dis- 
tance grid used by ACE to output distance-dependent quantities 
such as the size distribution was 100 values between 4 AU and 
400 AU. 

In many previous studies, the initial radial and size distri- 
butions of dust parent bodies — planetesimals — were taken 
in the form of power laws, with normalization factors and in- 
dices being parameters of the simulations. In this paper, we use 
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a different approach. To take advantage of our knowledge of the 
(largest) parent bodies, TNOs, we directly filled the (m, q, e)-bins 
at the beginning of each simulation with the objects of the "true" 
Kuiper belt. For comparison, we also made a run, where we pop- 
ulated the bins with known TNOs only (without debiasing). 

As already described, our "true" distribution contains only 
big objects with radii greater than ~ 5 km, whereas in reality 
small objects at all sizes down to dust must be present, too. 
If we started a simulation without smaller objects, the colli- 
sional cascade would take many gigayears to produce a note- 
worthy amount of dust and to reach collisional equilibrium. 
Accordingly, we have extrapolated the contents of the filled bins 
towards smaller sizes with a slope of ^ = 3.03 for objects be- 
tween 100m < i < 75km (in the gravity regime) and q - 3.66 
for objects smaller than 100 m (in the strength regime), follow- 
ing Fo'Brien & Greenberg (2003). Note that the adopted slope in 
the gravity regime is roughly consistent with Fig. [H The break 
at radi i of s everal t ens of kilometers reported bv Bernstein et al.l 
(l200l and lF^ase^ (l2009l) was not included. In the course of the 
collisional evolution, this artificial distribution corrects itself un- 
til it comes to a collisional quasi-steady state. The latter is as- 
sumed to have been reached, when the size distribution no longer 
changes its shape and just gradually moves down a s a whole as 
a resu lt of collisional depletion of parent bodies (iLohne et al.l 
I2OO8 '). We find that after £ 100 Myr a collisional quasi-steady 
state sets in for all solids in the strength regime (i.e., smaller 
than ~ 100 m). This particularly means that the system "does 
not remember" anymore the assumed initial distribution in this 
size range. 

We stress that the extrapolation of the observable EKB to- 
ward smaller sizes described above should not be misinterpreted 
as an attempt to describe the primordial size distribution of solids 
in the early EKB. The latter is set by the mechanism of the 
initial planetesimal accretion, which is as yet unknown. In the 
standard scenarios of "slow" accretion (e.g. lKenvon & BromlevI 
l2008i) a broad size distribution is expected. In contrast, "rapid" 
sce narios, such as the "primary accretion" mechanism proposed 
by ICuzzi et aL I (l200l or "graviturbulent" formation triggered 
by transient zones of high pressure ( Johansen et al. 2006) or 
by streaming instabilities (Joha nsen et al.i 2007) all imply that 
most of the mass of just-formed planetesimals was contained 
in 5 ~ 100 km bodies. Whatever mechanism was at work, and 
whatever size distribution the EKB in the early solar system had, 
in this study we are only interested in the present-day EKB. Thus 
the purpose of the extrapolation described above is merely to 
choose the initial size distribution across the sizes that would be 
as close to collisional equilibrium with the present-day EKB as 
possible. 
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3.2. Size distribution of dust 

FigurefTOldepicts the simulated size distribution of the EKB dust. 
We present three cases: for the debiased EKB without (top panel) 
and with P-R included (middle), as well as for the known EKB 
objects with P-R effect, for comparison (bottom). To explain the 
gross features of the size distributions shown in Fig. [TOl we 
introduce th e ratio o f radiation pressure to gravity, usually de- 
noted as j3 dBurns et al.lll979h . If a small dust grain is released 
after a collision from a nearly circular orbit, its eccentricity is 
e ^ /3(1 - py^ oc This implies higher eccentricities for 
smaller grain sizes. The orbits of sufficiently small grains with fi 
exceeding ~ 0.5 are unbound. Accordingly, the grain radius that 
corresponds toyS = 0.5 is commonly referred to as blowout limit. 
The blowout size for the assumed ideal material in the solar sys- 



Fig. 10. Size distribution of the Kuiper belt dust at different dis- 
tances. The vertical axis gives the cross section density per size 
decade. Top: the debiased EKB, without P-R effect; middle: the 
debiased EKB, with P-R effect; bottom: known EKB objects 
only, with P-R effect. 



tem is .Sbiow ~ 1 .2 ^m. Typically, the amount of blowout grains 
instantaneously present in the steady-state system is much less 
than the amount of slightly larger grains in loosely bound orbits 
around the star. This is because the dust production of the grains 
of adjacent sizes is comparable, but the lifetime of bound grains 
(due to collisions) is much longer than the lifetime of blowout 
grains (disk-crossing timescale). This explains a drop in the size 
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distribution around the blowout size which is seen in all three 
panels of Fig.fTOl 

Another generic feature of the size distribution is that it be- 
comes narrower at larger distances from the Sun. Were the par- 
ent bodies all confined to a narrow radial belt, the distribution far 
outside would appear as a narrow peak composed only of small, 
high-yS, barely bound grains sent by radiation pressure into ec- 
centric orbits with large apocentric distances. However, in the 
case of the EKB this effect is somewhat washed out, since the 
radial distribution of parent bodies themselves (mostly, of scat- 
tered objects) is extended radially, as discussed in Sect. 3.3. As 
a result, the size distribution even at relatively large distances 
(e.g. 76 AU) is a superposition of such a narrow distribution and 
a background broad distribution of particles produced at those 
distances directly. Only at largest distances, at which hardly any 
parent bodies are present (see 240 AU curve), the size distri- 
bution transforms to a predicted narrow peak adjacent to the 
blowout size. 

A direct comparison of our two simulations for the debiased 
EKB, without and with P-R effect, reveals some differences. One 
obvious — and expected — difference is the one between the 
24 AU curves. At this distance (and all the others inside the main 
belt) parent bodies are nearly absent; there are only some scat- 
tered TNOs, see Fig.[T] Accordingly, without P-R nearly no dust 
is present there. However, a substantial amount of small particles 
is present there in the P-R case, because these are transported 
there by the P-R drag. 

Outside ~ 30 AU, the size distributions without and with P- 
R show more similarities than dissimilarities. In particular, the 
maximum of the cross section at i = 2 \im is nearly the same. At 
sizes i ^ 1 mm, the curves roughly follow a classical Dohnanyi's 
law (cross section per size decade oc 5 "■^). The main difference 
is a dip of the size distribution in the region of the classical EKB 
that occurs at sizes of s - 100 [,im in the P-R case, which is 
easy to explain. The 100 \im grains in the classical EKB region 
stay in nearly-circular orbits, because their /3 ratio is small and 
radiation pressure-induced eccentricities are low. These grains 
are mainly destroyed in collisions with most abundant smaller 
grains, several i^im in size. In the non-P-R case, the latter grains 
have their pericenters within the classical belt. Thus the colli- 
sions are "grazing", the collisional velocities relatively low, and 
the collisional desctruction of 100 [xm grains relatively ineffi- 
cient. When the P-R effect is switched on, this changes. The P-R 
transport lowers the pericenters of smaller projectiles, and in the 
classical belt, they collide with 100 ^m grains at higher speeds, 
which enhances their destruction and produces the dip. Note that 
this effect is absent farther out from the Sun. At a 76 AU dis- 
tance, for example, the collisions between smaller grains and 
100 (,im particles always occur far from the pericenter. Thus the 
P-R effect has little influence on the collisional velocities, and 
the size distribution in the non-P-R and P-R cases is similar. 

However, the dust distributions computed with P-R effect, 
but for expected (i.e., debiased) EKB and known (not debi- 
ased) EKB, which are shown in middle and bottom panels of 
Fig-Hni exhibit a striking difference. The P-R effect has only a 
moderate influence on dust produced by the expected (debiased) 
EKB, but a strong one on dust generated solely by known TNOs. 
This needs to be explained. The debiasing procedure makes the 
EKB more densely populated, and the resulting increase in the 
dustiness shortens collision timescales to make them compara- 
ble with the P-R transport timescales. The resulting optical depth 
of the dust disk is such that it lies roughly between the collision- 
dominated and transport-dominated regimes. Without debiasing 
the parent body population, the dustiness of the disk is by two 



orders of magnitude lower, and so is the optical depth of the 
dust disk. At that optical depth level, the EKB dust disk would 
be transport-dominated below s :$ 100 ^m (but still collision- 
dominated at larger sizes). This is illustrated by the lowest panel 
in Fig. [To] that presents the size distribution of dust that would 
be produced by known TNOs. It is seen that the size distribu- 
tion in such a transport-dominated disk differs from that in a 
collision-dominated one qualitatively. From ~ 100 ^m down to 
blowout limit, the size distribution flattens and turns over. This is 
because the smaller the grains, the faster their inward P-R drift. 
This transport removes small grains from the collisionally ac- 
tive region and thus they are present in smaller amounts. As a 
result, the maximum of the cross section density shifts towards 
s ~ 100 [im particles. 

3.3. Radial distribution of dust 

Figure [TT] presents the radial distribution of dust parent bodies 
and their dust, the latter simulated without and with P-R trans- 
port. 

We start with a radial distribution of parent bodies, TNOs 
themselves, shown in Fig.fTTh. In contrast to Fig.|6l we plot here 
the total cross section of the TNOs instead of the mass they carry, 
because it is the cross section that characterizes the efficiency of 
TNOs as dust producers. Besides, we use the distance from the 
Sun instead of the semimajor axis as an argument. Specifically, 
we plot the cross section in the 80 km-sized TNOs, but the radial 
profile for larger objects look similar. As expected, the distri- 
bution peaks in the region of the main belt (40-50 AU), where 
about 90% of the cross section comes from the classical EKB 
objects. Outside ~ 60 AU, the cross section is solely due to scat- 
tered objects. The distribution of the latter is quite extended ra- 
dially, it is nearly flat over a wide distance range from x 35 AU 
to more than 100 AU. 

We now move to a discussion of the radial distribution of 
dust. As noted above, smaller grains with higher /3 ratios ac- 
quire higher orbital eccentricities. As the eccentricities of par- 
ticles slightly above the blowout limit are the highest, their ra- 
dial distribution is the broadest, whereas larger particles stay 
more confined to their birth regions. In Fig. [TTb this effect can 
be seen from how the curves gradually change from the largest 
(s = llO^im) to the smallest bound grains (2^im). The former 
essentially follow the distribution of the parent bodies, while 
the latter exhibit a more extended, flatter radial profile. Finally, 
blowout grains (e.g., those with s - 0.15 ^im have an oc r^^ dis- 
tribution, as expected for a set of hyperbolic orbits streaming 
outward from their birth locations. 

Including the P-R drag (Fig. [Tit ) does little with largest 
grains, but modifies the profile of smaller ones (~ 1 fxm) substan- 
tially. The P-R transport inward steepens their profile. Besides, 
small particles are now present in high amounts at smaller dis- 
tances, even where no parent bodies are present, in contrast to 
the case without P-R. 

Figure [TTB is the same as Fig. [TTb. but replotted in the same 
units as Fig. la in Moro-Martm & Malhotra (2003). Instead of 
treating the actual distribution of TNOs and their collisional 
dust production, they assumed a simplified narrow birth ring 
between 35 AU < a < 50 AU and a constant dust production 
rate. Nonetheless, the general behavior remains the same. At 
r a; 50 AU the number density begins to decrease rapidly but we 
still have some parent bodies outside 50 AU which produce dust, 
so that in our case the number density decreases more gently. 
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Fig. 11. Radial distribution of parent bodies ((a), from 64 km to 100 km in radius) and different-sized dust grains without (b) and 
with P-R drag (c). The vertical axis is the cross section density. The d panel is the same as c, but replotted in the same units as 
Fig. la in Moro-Mai-tm & Malhoti-a (.2003.) . 



3.4. Coupled size-radial distribution 

Another view of the EKB dust can be achieved by plotting its 
combined radial and size distribution (Fig.[T2ll. Besides present- 
ing the same salient features as those discussed in Sect. 3.2 and 
3.3, it emphasizes that radial and size distribution of dust in a 
debris disk are intrinsically coupled and cannot be treated inde- 
pendently of each other. 

3.5. Optical depth 

The radial profile of the normal geometrical optical depth is 
shown in Fig. [T3] In the case o f a disk dominated by the P-R 
effect, Strub be & Chian g (2006) calculated analytically the ex- 
ponent of the optical depth profile t oc r " to be a = 2.5 and 
or = in the outer and inner regions, respectively. Without P-R 
effect, i.e. for a collision-dominated disk, no dust is present inte- 
rior to the parent bodies, and the outer slope should be close to 
a = 1.5. 

These slopes are in qualitative agreement with our simula- 
tions (Fig. [T3T l. Taking the known EKB objects as dust sources 
and including the P-R effect, we find a nearly constant optical 
depth in the inner region and a slope of c x 3.0 in the outer 
disk, close enough to predictions for a transport-dominated disk. 
For the actual EKB dust disk and with the P-R effect taken into 



account, the result is intermediate between what is expected for 
transport-dominated and collision-dominated disks. This is seen 
from the inner profile which is gently decreasing inward (cf. 
Fig. 1 of IWvattll2005h . and from the outer slope of a a; 2.0. 
Finally, for the actual EKB dust disk, but with the P-R ef- 
fect switched off, the profile is the one expected for collision- 
dominated disks. The optical depth drops sharply inward from 
the main belt, whereas in the outer region the slope is or » 1.1. 
That it is somewhat flatter than the analytic value a - 1.5, traces 
back to a rather broad radial distribution of scattered TNOs that 
make a considerable constribution to the overall dust profile. 

For the debiased EKB and with P-R effect included, the nor- 
mal optical depth peaks at * 40 AU at a level of ^ 6 x 10"^. 
Besides the normal optical depth shown in Fig. [13] we have cal- 
culated the in-plane optical depth. Our result, t < 2 x 10"^ out- 
si de 30 AU, is at least by a factor of 4 higher than the estimate 
bv lSteriil ([19961) who found r = 3 x 10"^ . . . 5 x 10 ^'. 

We finally note that the dust production rate for the expected 
EKB in collisional equilibrium was calculated to< 1.7x10^ gs"' 
and is three times hi gher th an the predicted rate of 5 x 10^ gs"' 
from lLandgraf et al.1 (l2002h on the base of in situ measurements 
of Pioneer 10 and 1 1 . 
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Fig. 12. The distribution of the cross section density of the EKB 
dust as a function of distance and grain sizes. The panels are as 
in Fig. [m 



4. Spectral energy distribution 

4.1. Parameters and materials 

The equilibrium temperatures of dust and their thermal emission 
w ere calculated in a standard way as described, for instance, 
in iKrivov etal ] (I2008h . In these calculations we computed the 
solar photosphe ric spectrum with the NextGen grid of models 
dHauschildt et al. 1999), assuming a G2V star of solar metallic- 
ity. To get a rough idea of how the thermal flux is affected by 



10- 



^ 10^ 



10 ' 



10" 



w/o PR, expected 
with PR, expected 
with PR, l^nown 



a = 1.1 

•,a = 3.0 ^^a = 2.0'' 



,,a= 1.5 



a = 2.S 



20 



30 40 50 60 7080 100 
r[AU] 



200 



400 



Fig. 13. Normal optical depth of the Kuiperbelt with and without 
P-R drag. Uranus and Neptune are shown for orientation, but 
were not included in the simulations. The optical depth for the 
known EKB is amplified by a factor of 100. 



(unknown) chemical composition of the dust, we tested three 
different cases: ideal m aterial (blackbody abso rption and emis- 
sion), astrosilicate from iLaor & Drain3 (Il993h . and "dirty ice". 
The latter is defined as a mixtu r e cons isting of a water ice matrix 
(with constants from IWarrenI (119841) ') c ontaminated with 10% 
volume fraction of astronomical silicate (Laor &Draine"1993|). 
This is similar to what, for instance, Yamamoto & Mukai ( 199q) 
adopted in their calculation of thermal emission of the EKB dust. 
The refractive indices of the adopted ice-silicate mixture were 
calculated by means of the Maxwell-Garnett theory. This also 
explains our using 10% silicate fraction: 10% is an upper limit 
for which the effective medium theory still provides accurate re- 
sults (Koloko lova & GustafsonI 120011) . With a standard Mie al- 
gorithm, we then computed the absorption efficiency Qabs as a 
function of size and wavelength. 

4.2. Calculated flux 

The resulting SEDs of the simulated EKB dust disk, as it would 
be seen from a 10 pc distance, are presented in Fig. [141 Three 
different curves correspond to three materials described above. 
Pure astrosilicate and the ice-silicate mixture produce SEDs of 
similar shape and height, peaking at 50-70 ^m with a maximum 
flux at a level of several mJy. On the "dirty ice" curve, a typi- 
cal water ice feature at ~ 60 (xm is seen. This feature may have 
been detected in th e disk of a young debris disk star HD 181327 
(IChen et alJl2"008h . The fact that the feature is located near the 
SED maximum may help finding water ice in other debris disks 
by future observations. A comparison with the blackbody curve 
readily shows that it departs from the two others significantly. 
The maximum moves to A ~ 200 (xm and the (sub)-mm flux 
becomes by two orders of magnitude higher than in the two 
other cases. This result confirms earlier conc lusion (see, e.g. 
'Yamam oto & Mukailll998l; iKrivov et al.]|2008l) that the black- 
body assumption is probably too crude and should not be used 
in modeling the thermal emission of debris disks. 

It is important to compare the calculated flux of the EKB 
dust with measurements of other known debris disks. For this 
purpose, we looked at photometric data for other Sun-like stars. 
We have chosen three stars which possess well-known, bright 
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Fig. 14. Infrared excess of the Kuiper belt dust, (a): the debiased 
EKB, without P-R effect; (b): the debiased EKB, with P-R ef- 
fect; (c): known EKB objects only, with P-R effect (amplified by 
a factor of 100); (d): the ratio of (b) and (c). Three curves in each 
panel (a) to (c) are based on the same ACE runs, but the ther- 
mal emission was calculated for three different materials: black- 
body, astrosilicate, and "dirty ice". For comparison, observed 
SEDs of three other old Sun-like stars are shown with symbols. 



Photometric data for these stars are adapted from Greaves et al 
(2004); Sheret et al. (2004); Beichman et al. (2006); Chen et^ 
(2006); Moor et al.1 (120061) : ILawler et all (120091) : ITanner et al 
(.2009) . 



G8V, 4 pc, 7Gyr). Their SEDs, normalized to r = 10 pc, are 
overplottted in Fig. [14] A comparison shows that their SEDs are 
similar in shape to the simulated SED of the EKB dust disk (for 
one of the two realistic materials, not blackbody). The maximum 
of the SEDs of the selected debris disk stars lies at ~ 100 \im, i.e. 
at slightly longer wavelenghts than the maximum of the EKB 
dust flux. This (moderate) difference in the peak wavelengths of 
different SEDs may be caused by a choice of chemical composi- 
tion of the disks, by different extentions of EKB and other debris 
disks, or both. 

The major difference between the EKB dust disk and extra- 
solar debris disks is, of course, the absolute level of the thermal 
emission fluxes. For dust maintained by the debiased EKB (pan- 
els a and b), the fluxes from far-IR to sub-mm are about one to 
two orders of magnitude lower that those of the reference stars. 
Obviously, this traces back to a much lower total mass of the 
EKB compared to that of the extrasolar Kuiper belts. The rea- 
sons f or this d ifference are currently a matter of debate. For in- 
stance, iBooth et al. (2009) argue that typical Kuiper belts around 
stars with observed debris disks may not have undergone a major 
depletion phase due to a "Late Heavy Bombardment". Currently, 
it is not known whether low-mass disks at the EKB level exist 
around other stars at all and, if they do, how common they are. 
This question can only be answered with future, more sensitive 
observations (for instance, with Herschel Space Observatory). 

It is interesting to compare fluxes from the expected EKB 
dust disk (Fig. [T4b ) with those from the dust disk that would 
stem from the known TNOs, without debiasing. The correspond- 
ing SEDs, depicted in Fig. [T4b . are completely different. Apart 
from an expected reduction of the fluxes by almost two orders 
of magnitude, two other effects are seen. One is that the as- 
trosilicate and "dirty ice" SEDs come much closer to each other 
and to the blackbody curve. As a result, their maxima shift to 
longer wavelengths of ~ 100|,im. Another effect is a disapper- 
ance of the water ice feature. Both effects are easily explained 
by the major differences between the size distribution of dust 
of the known and expected EKBs (cf. middle and bottom pan- 
els in Fig. [Tot . In the dust disk of the known EKB, the cross 
section, and thus the thermal emission, are dominated by grains 
~ 100 nm in size. Such big grains behave as black bodies and 
do not produce any distinctive spectral features. We also present 
the ratio of fluxes from the dust disks of the expected and known 
EKBs (Fig.[T4}l). The aforementioned reduction of the fluxes is 
seen here even better. The reduction factor, ~ 10 at wavelengths 
longward of ~ 100 |,im, increases further at shorter wavelengths 
for ice-silicate and especially astrosilicate grains. Obviously, this 
is also a consequence of the change in the size distribution dis- 
cussed above. 

For three selected wavelengths of 70 ^m, 100 ^m, and 
160 [,im, the fluxes calculated for the "known" and "expected" 
EKBs and for all three materials are also tabulated in Table [5] 
Solar photosphereic fluxes are given for comparison. These val- 
ues are used below to assess the detectability of debris disks with 
Herschel. 



4.3. Detectability of "Kuiper beits" 



excesses (measured at many wavelengths from mid-IR to sub- 
mm), which have spectral classes not too far from solar G2, 
and which, like the Sun, are rather old (^ 1 Gyr). These are HD 
48682 (i/r^ Aur, GOV, at 17 pc distance, estimated age 0.6-9 Gyr), 
HD 10647 (qi Eri, F8V, 17 pc, ~ 2 Gyr), and HD 10700 (r Cet, 



We would like to estimate the level, down to which Herschel can 
detect faint debris disks, and then to compare that level with our 
EKB dust disk models. 

The sensitivity limit of the PACS instrument of the Herschel 
telescope at 100 ^im in the scan-map mode for one hour expo- 
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Table 5. Calculated flux of the known and expected Kuiper belt 
dust (including the P-R-effect) and photospheric flux of the Sun 
at different wavelengths from a lOpc distance. The fluxes are 
given in mJy. 



sure time is 2.4 mJy at a 5cr uncertainty leveQ, or 0.48 mJy at 
la. The background noise, of course, is highly variable from 
one star to another. For 133 target stars of the Open Time 
Key Program DUNES (DUst around NEarby Stars), its aver- 
age value at 100 (im is 0.53 mJy at Icr. Combining the instru- 
ment and the background noise leads to a limiting flux as low as 
5 VO.48^ + 0.53^ = 3.6 mJy for one hour exposure time at a 5cr 
uncertainty level. Assuming that the 100 (,im flux is proportional 
to the total mass of a debris disk, we can conclude that an extra- 
solar EKB analog with a mass M > 1 .7Mekb will be detectable 
from a distance of lOpc, assuming a "dirty ice" as a dust com- 
position. For silicate dust, the detectability limit would go down 
to 0.9Mekb- Additional uncertanties in the stellar photospheric 
flux (x 2% of 22 mJy, or 0.4 mJy) increase these values slightly 
to s; 1.9Mekb and x I.OMekb, respectively. 

5. Summary and conclusions 

In this paper, we attempted to construct a model of the solar sys- 
tem debris disk — Kuiper belt objects and their collisional de- 
bris. This was done in two major steps. First, we developed a new 
algorithm to remove the two largest selection effects from TNO 
observations, the inclination and the distance biases. Applying 
it to the database of 1260 known TNOs, we derived new esti- 
mates for the global parameters of the "true" Kuiper belt: its 
mass, orbital and size distribution. Second, treating the debiased 
populations of EKB objects as dust parent bodies, we employed 
our collisional code to simulate their dust disk. This has resulted 
in estimates of the size and radial distribution of dust in the 
transneptunian region, as well the optical depth of the Kuiper 
belt dust disk. We have also calculated the thermal emission of 
the Kuiper belt dust and compared the expected infrared fluxes 
with those of known debris disks around other Sun-like stars. 

We draw the following conclusions: 

1 . The total mass of the EKB, including classical, resonant, and 
scattered objects, amounts to 0.12Me. Without scattered ob- 
jects, the mass reduces to 0.05 M^. These values are by about 
one order of magnitude higher than the mass of the known 
TNOs. 

2. The dust disks that would be produced by the known EKB 
and by the expected (i.e., debiased) one are shown to have 
distinctly different properties, as detailed below. Much of the 
difference in the size and radial distribution of dust of the two 
EKBs traces back to the increased total mass of TNOs after 
the debiasing. 

3. The slope of the differential size distribution in the EKB in 
the size range 50 km < s < 170 km (absolute magnitudes 

^ http://herschel. esac. esa. int/Docs/PACS/html/ch03s05. htmWsec- 
photo-sensitivity 



8.9 > H > 6) is found to be ^ = 4.3 ± 0.2 for the classical 
objects, 5.1 ±0.1 for resonant ones, and 2.8±0.1 for scattered 
TNOs. This results in 4.4 ± 0.2 for classical and all resonant 
objects together. For the entire EKB (classical, resonant, and 
scattered TNOs), we find 3.6 ± 0.1. 

4. If the debiased populations of TNOs are taken as dust 
sources, the dustiness of the disk is high enough to make col- 
lisional timescales comparable to, or even shorter than, the 
radial transport timescales. Thus the predicted Kuiper belt 
dust disk rather falls to the category of collision-dominated 
debris disks, to which all debris disks detected so far around 
other stars belong, albeit it is already close to the "bound- 
ary" between collision-dominated and transport-dominated 
disks. The simulated size distribution in the EKB dust disk 
shows a sharp cutoff at the radiation pressure blowout radius 
of ~ 1 ^im. The cross section dominating radius is several 
times larger. 

5. For comparison, we have also considered the dust disk that 
would be produced solely by known TNOs as parent bodies. 
That disk would fall into the transport-dominated regime, 
where the Poynting-Robertson drift timescales are shorter 
than the collisional timescales (at dust sizes). The size dis- 
tribution in such a disk is dramatically different from the one 
expected in collision-dominated disks. While a cutoff at 1 |,im 
remains, the distribution between 1 \im and several hundreds 
of |,im is nearly flat, and the cross section is now dominated 
by much larger grains, ~ 100 ^m in radius. 

6. The radial distribution of the TNOs themselves peaks in the 
region of the main belt (40-50 AU), where about 90% of the 
cross section comes from the classical EKB objects. Outside 
~ 60 AU, the cross section is solely due to scattered objects. 
The distribution of the latter is quite extended radially, being 
nearly flat over a wide distance range from » 35 AU to more 
than 100 AU. 

7. The radial distribution of dust grains with radii ^ 10(,im 
is similar to the distribution of the parent bodies described 
above. At smaller sizes, the radial profile gets progressively 
broader with decreasing radius, which is a classical effect of 
radiation pressure. 

8. The maximum normal geometrical optical depth is reached 
at the inner edge of the classical belt, ^ 40 AU, and is es- 
timated to be » 6 X 10"''. Outside that distance, it falls 
off as r^^. This slope is roughly intermediate between the 
slope predicted analytically for collision-dominated (r"'^) 
and transport-dominated (r^^^) disks. An upper limit of the 
in-plane optical depth is set to a; 2 x 10"^ outside 30 AU. 

9. For comparison, the normal optical depth of the dust disk 
produced only by known TNOs would fall off outside the 
classical EKB as r oc r"^, and the in-plane optical depth 
would be by two orders of magnitude lower. 

10. The estimated thermal emission flux from the EKB dust disk 
that would be observed from a lOpc distance strongly de- 
pends on the assumed dust composition. For two "realistic" 
materials probed, astrosilicate and ice-silicate mixture, the 
SED would peak at 40-60 (,im. The maximum value is at a 
level of several mJy (at the solar photospheric flux at the 
same wavelengths is several tens of mJy), which is about 
two orders of magnitude smaller than for the brightest known 
debris disks around other solar-type stars, normalized to the 
same distance. For a non-debiased EKB, the fluxes would be 
another two orders of magnitude lower. 

1 1 . With the Herschel Space Observatory, it may be possible to 
detect faint debris disks around other stars, which are nearly 
as faint as the EKB dust disk. For observations with the 
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PACS instrument at 100 [^im in a scan-map mode with one 
hour exposure, assuming an average background noise, the 
minimum mass of a Kuiper belt analog detectable at 5cr from 
lOpc is estimated as a: 1 . . . 2Mekb- 

After the submission of this pap er, we became aware of 
another work (lKuchner&Starkll20Ta ) that addresses the same 
problem of predicting the EKB dust disk properties and, like 
this paper, employs a catalog of known KBOs to model the 
dust production. Instead of collisional grinding simulations, 
iKuchner & StarkI make use of N-body integrations supple- 
mented with a "collisional grooming" algorithm. Despite the dif- 
ference in the simulation methods, their results for the radial dis- 
tribution of the EKB dust are in a very good agreement with ours 
(cf. our Fig, nn and their Figs. 4-5). The same holds for size dis- 
tributions at lower optical depth levels where they find the same 
turn-off in the size distribution as we do (cf. our Fig. \W\ and 
their Fig. 6). As a result, the size of cross section-dominating 
grains in the densest parts of disks becomes larger. Since these 
larger grains are cooler, the maximum of the SED should move 
to longer wavelengths with decreasing r (cf. our Fig.[T4]and their 
Fig. 9). 



6. Discussion 

The main result of this paper is a set of new estimates of the 
parameters of the "true" (debiased) Kuiper belt and its debris 
dust. As all models, ours involves uncertainties. Here we discuss 
three major ones. 

7. Parent bodies. Although we believe that our debiasing al- 
gorithm is more accurate that the ones used before, it is possible 
that the results over- or underestimate, for instance, the mass of 
the EKB by a factor of a few. The main reason is an incomplete 
list of surveys used for debiasing. In particular, observations that 
covered higher ecliptic latitudes, omitted in Tab. [T] could alter 
the results. We have identified 90 objects in the MPC (excluding 
Centaurs) that have been detected beyond 10° ecliptic latitude. 
Of these, 29 objects have been found beyond 20°. To check how 
non-including the surveys that covered higher latitudes may af- 
fect the results, we have made additional tests. We have re-run 
our code several times, each time artificially changing s in three 
randomly chosen surveys of Tab. [1] to 20°. This changed the de- 
tection probability of single objects typically by a few tenths of 
percent only. The change in the inferred total EKB mass turned 
out to be larger, but still minor, 4%. 

Another difficulty is related to the scattered objects with 
their large eccentricities and inclinations. Their average detec- 
tion probability is very low, the resulting "mass amplification 
factor" very high (» 25), so that the estimated masses of objects 
can considerably diff'er from the "true" ones. 

The same applies to the treatment of the smallest ob- 
jects discovered so far, a few kilometers is radius. If, for in- 
sta nce, the break t he siz e di stribution at s ^ 30 km inferred 
bv lBernstein et al.1 (|2004 and lFrased (l2009l) were included, this 
would reduce the overall amount of dust. This would make the 
EKB dust disk "more transport-dominated" with all the conse- 
quences disussed above. 

Being aware of such caveats, throughout the paper we always 
compare the known (not debiased) EKB with the expected (de- 
biased) EKB. This applies to the EKB itself (Sect. 2), as well as 
to its dust (Sect. 3-4). We expect that the results obtained here 
for the known and debiased EKB would at least "bracket" the 
reality. 



2. Material composition of solids. Realistic simulations 
would necessitate good knowledge of tens of parameters, such 
as the bulk density, porosity, shape, tensile strength, optical con- 
stants of solids in the EKB, and all this over an extremely broad 
range of sizes: from hundreds of kilometers (large TNOs) down 
to a fraction of a micrometer (tiniest dust grains). In this paper, 
we made a large set of simplifying assumptions, ignoring, in par- 
ticular, any dependence of all these parameters on size, although 
such dependencies are expected. We only probed several diff'er- 
ent materials when calculating the thermal emission of dust in 
Sect. 4, and saw how sensitive the results are to the adopted ab- 
sorption and emission efficiencies of grains. There are no reasons 
to think that variation of other parameters, for instance mechan- 
ical properties of diff'erent-sized solids in the collisional model- 
ing, would alter the results to a lesser extent. 

3. Modification of dust disk by Neptune. As mentioned 
above, effects of giant planets on the dust distribution, of 
which those by Neptune are the largest, were not included 
in our simulations. These effects are diverse and can alter 
the results substantially. Neptune is expected to capture some 
grains stemming from non-resonant KBOs and d rifting inward 
by P-R effect into mean-motion resonances (e.g. iLiou & Zookl 
119991; iMoro-Martm & Malhotral |2002 |). Similarly, sufficiently 
large grains that derive from resonant TNOs must stay locked in 
such resonances from the very beginning for considerable time 
periods (e.g. Wyatt 2006; Krivov^elal. 2007). This will cause an 
enhancement of dust density at resonance locations, which was 
not modeled here. Furthermore, the distribution of dust will ex- 
hibit azimuthal clumps instead of being rotationally-symmetric, 
as was assumed in our ACE simulations. Next, many of the 
grains are in Neptune-crossing orbits, or reach such orbits in the 
process of their inward drift. The probability of close encounters 
with Neptune is not negligible, especially for larger grains that 
drift inward more slowly. Encounters would result in partial trun- 
cation of the dust disk at the planet's orbit, i.e. in a drop of dust 
density at 30 AU. Some of the grains will be ejected outward, 
which may affect the collisional balance of the whole disk. 

For these reasons, the simulation results presented here 
should be treated as preliminary. We can hope that future, deeper, 
TNO surveys would lead to a better knowledge of the global pa- 
rameters of the EKB acting as dust parent bodies. Further mod- 
eling work, perhaps with other methods, should help quantify 
the effects of Neptune on the dust distribution. Choosing the 
most adequate mechanical and optical properties of KBOs and 
their dust is probably the main difficulty. Eventually, many of 
the uncertainties could be minimized, and the models can be 
verified, if dust in the transneptunian region is detected in situ 
by dust detectors ab oard space missions, such as New Horizons 
(iHoranvi et alJl20()8 ). On any account, we consider this study as 
a reasonable starting point in developing more accurate models 
of the debris disk of our own solar system. 
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